Radiative decay of the X(3872) as a mixed molecule-charmonium state in QCD 

Sum Rules 
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We use QCD sum rules to calculate the width of the radiative decay of the meson X(3872), 
assumed to be a mixture between charmonium and exotic molecular [cqJIqc] states with J^'~^ — l^"*". 
We find that in a small ranee for the values of the mixing angle, 5" < e < 13^ we get the branching 
ratio T{X — >■ J/il)'y) /T{X — >■ J/tpTT'^TT~) = 0.19±0.13, which is in agreement, with the experimental 
value. This result is compatible with the analysis of the mass and decay width of the mode ,J/ip(mi) 
performed in the same approach. 
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I. INTRODUCTION 



The X(3872) state has been first observed by the 
Belle coUaboration in the decay B+ X(3872)i4:+ 
J/'f/'TT+TT^if + Jl|, and was later confirmed by CDF, DO 
and BaBar ,31 ■ The current world average mass is 
mx = (3871.4±0.6) MeV, and the width is F < 2.3 MeV 
at 90% confidence level. Babar collaborations reported 
the radiative decay mode X(3872) — )■ 7J/V' [1, which 
determines C = +. Belle Collaboration reported the 
branching ratio: 



T{X J/iPt:+it-) 



= 0.14 ±0.05. 



(1) 



Further studies from Belle and CDF that combine an- 
gular information and kinematic properties of the tt+tt^ 
pair, strongly favors the quantum numbers J^*^ = 1++ 
or 2 0, Between these quantum numbers, a re- 

cent BaBar measurement favors the J^'-^ = 2 ^ assign- 
ment 0- However, established properties of the X(3872) 
are in conflict with this assignment and, there- 

fore, in this work we assume the quantum numbers of 
the X(3872) to be J 



PC _ 2^-1-+ 



The interest in this new state has been increasing, since 
the mass of the X(3872) could not be related to any char- 
monium state with the quantum numbers J^*^ — 1++ in 
the constituent quark models [lo| . indicating that the 
conventional quark-antiquark structure should by aban- 
doned in this case. Another interesting experimental 
finding is the fact that the decay rates of the processes 
X(3872) J/V'TT+Tr-TrO and X(3872) J/V'tt+tt- are 
comparable [3j: 



T{X J/TAvr+Tr-vr") 
T{X -^J/iPtt+t:-) 



1.0 ±0.4 ±0.3. 



(2) 



This ratio indicates a strong isospin and G parity vi- 
olation, which is incompatible with a cc structure for 
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X(3872). The decay X — > J/il>uj was also observed by 
BaBar Collaboration j7] at a rate: 



B{X J/iZ-tt+tt-ttO) 

B{X ^ J/V'TT+TT-) 



= 0.8 ±0.3, 



(3) 



which is consistent with the result in Eq. 

The isospin violation problem can be easily avoided 
in a multiquark approach. In this context the molecu- 
lar picture has gained attention. The observation of the 
above mentioned decays, plus the coincidence between 
the X mass and the D*°D° threshold: M{D*°D°) = 



(3871.81 ±0.36) MeV [ll|, inspired the proposal that the 
X(3872) could be a molecular {D*°D° - D*°D°) bound 
state with small binding energy [H The 
molecule is not an isospin eigenstate and the rate in 
Eq. ^ could be explained in a very natural way in this 
model. 

Although the molecular picture is gaining attention 
with studies indicating that it can be a suitable descrip- 
tion for the X(3872) structure p^, there are also some 
experimental data that seem to indicate the existence of 
a cc component in its structure. In ref. [l5l | , a simulation 
for the production of a bound D^D*^ state with biding 
energy as small as 0.25 MeV, reported a production cross 
section that is an order of magnitude smaller than the 
cross section obtained from the CDF data. A similar re- 
sult was obtained in ref. [l6| in a more phenomenological 
analysis. However, as pointed out in ref. [T^, a consis- 
tent analysis of the D'^D*^ molecule production requires 
taking into account the effect of final state interactions 
of the D and D* mesons. 

Besides this debate, the recent observation, reported 
by BaBar [18], of the decay X(3872) -> '0(2S')7 at a 
rate: 



B{X^4,{2S)-i) 
B{X ^ ^7) 



= 3.4 ± 1.4, 



is much bigger than the molecular prediction [IJ 

r(x^ V(25)7) 



V{X ^ ^7) 



4 X 10 



-3 



(4) 



(5) 



Another interesting interpretation for the -X_(3872) is 
that it could be a compact tetraquark state [20l - |23| . In 
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particular, Terasaki [22] argues that with a tetraquark 
interpretation the ratio in Eq. ([T]) could be easily ex- 
plained. 

In Ref . j23| the QCDSR approach was used to study 
the X structure including the possibility of the mixing 
between two and four-quark states. This was imple- 
mented following the prescription suggested in (25l | for 
the light sector. The mixing is done at the level of the 
currents and is extended to the charm sector. In a dif- 
ferent context (not in QCDSR), a similar mixing was 
suggested already some time ago by Suzuki [l^. Phys- 
ically, this corresponds to a fluctuation of the cc state 
where a gluon is emitted and subsequently splits into a 
light quark-antiquark pair, which lives for some time and 
behaves like a molecule-like state. The possibility that 
the X(3872) is the mixing of two-quarks and molecu- 
lar states was also considered to investigate the radiative 
decay in the effective Lagrangian approach [2(?|, and to 
explain the data from BaBar and Belle 27] using a 
Flatte analysis [H. 

In this work we will focus on the radiative decay 
X(3872) — > J/il'l- We use the mixed two-quark and four- 
quark prescription of Ref. [2^ to perform a QCD sum rule 
analysis of the radiative decay X(3872) — )■ J/V'7. 



II. THE MIXED TWO-QUARK / FOUR QUARK 
OPERATOR 

The mixed charmoniuni-molecular current proposed 
in Ref.jl^l will be used to study radiative decay of the 
X(3872) in the QCD sum rules framework. 

For the charmonium part we use the conventional axial 
current: 

= Ca{xhf,-f5Ca{x). (6) 

The D D* molecule is interpolated by [H-Illl: 
1 



^ ^ {qa{x)-f5Ca{x)cb{x)-i^,qb{x)) 
{qa{x)jf^Caix)cb{x)-f5qb{x)) 



(7) 



As in Ref. |25| we define the normalized two-quark cur- 
rent as 

1 , 



,■(2?) 



6V2 



(8) 



and from these two currents we build the following mixed 
charmonium-molecular current for the X(3872): 

J^(x) = sin(0)j(4«) {x) + cos(0)j(2^) (x). (9) 

Following Ref. f23] we will consider a D^^D*^ molecular 
state with a small admixture of D'^D*~ and D^D*^ 
components: 



jnix) = COS a JJ;^ (x) + sin aJ^{x), 



(10) 



with J^{x), {q — u,d), given by the mixed two- 
quark/four-quark current in Eq. 



III. THE THREE POINT CORRELATOR 

In this section we use QCD sum rules to study the ver- 
tex associated to the decay X (3872) J/4'1- The QCD 
sum rules approach |32 - l34l | is based on the principle of 
duality. It consists in the assumption that a correlation 
function may be described at both quark and hadron 
levels. At the hadronic level (the phenomenological side) 
the correlation function is calculated introducing hadron 
characteristics such as masses and coupling constants. 
At the quark level, the correlation function is written in 
terms of quark and gluon fields and a Wilson's opera- 
tor product expansion (OPE) is used to deal with the 
complex structure of the QCD vacuum. 

The QCD sum rule calculation for the vertex 
^(3872) J 14'! is centered around the three-point func- 
tion given by 



^^iua{p,p' ,q) 



d^xd^ye'P'-^ e'''-yU^,^{x,y), (11) 



with 



Tl,.aix,y) = {0\T[jfix)j2iy)j^\0)]\0), (12) 

where p — p' + q and the interpolating fields are given 
by: 



^ eqqji^q, 

q—u,d,c 



(13) 



(14) 



with Gq — for quarks u and c, and = — for quark 
d (e is the modulus of the electron charge). The current 
is given by the mixed charmonium-molecule current 



in Eq. ([T0|. 

In our analysis, we consider the quarks u and d to be 
degenerate, i.e., m„ = and (uu) = (dd) = {qq)^ then 
by inserting the mixed current (jlOp in Eq. (|12p , we arrive 
at the following relation for the correlator 



e sm ( 



n^i/a {x,y)^ ^^^^ (2 cos a - sin a)n™°^ (a;, y) 



3 

e(gg) 

6^2 



cos0(cos a -f sinQ;)n^'j,„(a;, y) . 



(15) 



The relation for the correlator is written in terms of the 
charmonium and molecule contributions. For the char- 
monium term we have 

K^Ax.y) = (0|r[j;C(x)j7(y)jl(2)^(0)]|0), (16) 

and the molecular term is given by 

K:l{x,y) = {mjtix)j2{y)j^y\Q)m, (17) 

with ja"^^ and j'"a'^^ given by Eqs. (jS]) and (O respectively. 
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We now proceed to the calculation of both charmo- 
nium and molecular contributions in in the OPE side. 
By inserting the currents of the two-quark component, 
J/ipi and photon, respectively defined in Eqs. 
and (|14p . in Eq. (jl6p . we obtain for the charmonium con- 
tribution the following relation: 

2 

n^ta(a^:y) = -3 Tr [l^Sl^[x - y)j^S^^{y)'yal5S^a{-x) 
+ lt.SlAx)lc.l5Sl^{-y)l.Sl^{-x + y)] , (18) 



where Slf^{x - y) = {0\T[qa{x)qb{y)]\Q) is the full propa- 
gator of the quark q (here a, 6, c are color indices). 

For the molecular contribution we use the four-quark 
current defined in Eq. ([7]), as well as the currents for the 
J/ij} and the photon. Inserting these currents in Eq. ([T7|. 
we get 



ttOPE 



{x,y) 



V2 



Tr 



lf^S2'aixh5Sl^,i-y)-f:, X 



xSl,{y)^o.SL,{-x) - l^.Sl.MlA,{~yh. x 



-xSl,,{y)^,SUk-p') 



(19) 



To evaluate the phcnomenological side of the sum rule 
we insert, in Ec dTTl) . intermediate states for X and J/il). 
We use the following definitions: 

mjtW))^m^u^M)-^ (20) 

{X{p)\j^\Q) = (cosa + sina)A,e;(p), (21) 

where the meson-current coupling parameter is extracted 
from the two-point function, and its value was obtained 
in Ref. [13: A, = (3.6 ± 0.9) x lO^^ QeV^ We obtain 
the following expression: 

(cos a -I- sina)Aqm^/^e^(p')e* (p) 



n^r(p,p',9) = -- 



X {i^{p')\j2\X{p)) 



(22) 



The remaining matrix element can be related to the one 
that describes the decay X ^ ^J/ip: 



mp'mmxip)) = ieZ{q)M{x{p) ^ i{q)Jiw)) , 

(23) 

and we can define [2^ 

M{X{p) ^ i{q)J/i;{p')) = es-^P^e^ApK^P'mq) x 
x-^{Ag^xgaKP ■ q + Bg^xPKQa + Cga^pxqp,), (24) 

where A,_B,C are dimensionless couplings. Using this 
relation in Eq. (j22p . we can write the phcnomenological 



side of the sum rule as: 



ze(cosQ! + sin a) Xqiriti,/^ 
m\{p'^ - ■m\){p''^ - m^) 



X e"'^''"<7.p • qA -f e^-^-'^p'xqaqo.B - q^q^p'^C 



p ■ q 



e^'^^y'xqAqc+p'JiA + B) 



p ■ q 

2 

X 



(25) 



In the OPE side we work in leading order in as and 
we consider condensates up to dimension five, as shown 
in Fig. [TJ In the phenomenological side, as we can see 
in Eq. ([25|) . there are five independent structures. We 
choose one convenient structure to determine each one 
of the couplings A,B,C in Eq. ([M)) . Taking the fimit 
p^ — p'^ = —P^ and doing a single Borel transform to 
Af^, we arrive at a general formula for the sum 
rule for each structure i: 



G 



n 



(OPE) 



(M2,Q2 



(26) 



where = —q^ and Hi{Q^) give s the contribution of 
the pole-continuum transitions |35l437j | . In the following, 
we show the expression of the sum rules for the three 
structures that we have chosen to work. 



1. Structure 1: e"'""'g<. 

The RHS of the sum rule for the structure e°'^'^"'qa 
(structure 1) have both charmonium and molecule con- 
tributions: 

sin 0(2 cos a — sin a) 



nr^(M^g^) = -(g-g) 



3g4 



xnf«(Af2, g2) + ^(cosa -I- sina)nf (M^, 



2g2 



where the molecular contribution is given by 



(27) 



3g2 



X \ 1 



4m? / 1 









du e 


J im^ 




'1 _^ 






u I 



l + 3a 



2 2 rl 



g o,(l-a)M-' 



(28) 



and the charmonium contribution is 

jccr,,2 - ds due-^ 



nf(A/2,g2 



Va 



(29) 
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(f) 



(g) 



FIG. 1. Diagrams which contribute to the OPE side of the sum rule. Diagrams (a) to (e) contribute to the molecule term of 
the OPE; (f) and (g) contribute to the charmonium term. 



where A = A(s, t, u) — + + — 2st — 2su — 2tu, and 
t = -Q2 ^ 0. 

In the above expressions the parameters Sq = {mx + 
As)^ and uq — {m^ + A^)^ are the continuum thresholds 
for X and J/?A respectively. The limits of the integral in 
u are: 



u± 



(30) 

The integrals in s and u also obey the following condi- 
tions: 



t < u, 4m^ < So . 



(31) 



Since the photon is off-shell in the vertex XJ/ip'^ it is 
required the introduction of form factors. Then in the 
left hand side of the sum rule, we define the function 
Gi{Q^), which is related to the form factor A{Q'^) as: 



3\/27r^ (cos a + sin a)Xqm^f^ 



A[Q^). (32) 



2. Structure 2: e'^'^^p'^p'^qx 



The RHS of the sum rule for the structure e^'"^^p'^p'^qx 
(structure 2) has only molecular contribution: 



nOPE(Af2,g2) 



da 
a 



a - 



e -(i-=)Af-(33) 



In the left hand side of the sum rule we define the 
function G2(Q^), which is related to the sum of form 



factor A{Q^) + B{Q^) as 
G2{Q 



2, _ 2,^2-^ V2tt^ {cos a + Sin a)\qm^f^{A{Q^) + B{Q^)) 



sin 9{2 cos a — sin a)mj^ (i^x ~ ™^ ) 



(34) 



3. Structure 3: e^^^^p', 



Pxq^QiJ. 



The RHS of the sum rule for the structure e"'^-*"^-"' 



(structure 3) has only charmonium contribution: 



^op^=(A^^g2) = (g<z) / ds 

J Ami 



dU e TirTTT X 



A3/2 



tu + TO^(— s + t + u) + 



3 s t u(— s + t + u) 



A 



(35) 



The integrals in this equation obey the same relations 
and conditions defined for the Eq. (|29p . 

In the left hand side of the sum rule we define the func- 
tion G3{Q^), which is related to the form factor C{Q^) 
as: 



6V2'K'^Xqm^f^ 
cos (m^ — m^) 



cm- 



(36) 



IV. NUMERICAL ANALYSIS 

The sum rules are analysed numerically using the fol- 
lowing values for quark masses and QCD condensates 
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(c) 

FIG. 2. Values of the functions obtained by varying both and M'^: a) Gi(Q^), b) G2{Q^) and c) GsCQ^). 



[23|,|38j, and for meson masses e decay constants: 

mdrric) = (1.23 ±0.05) GeV, 
(qq) = -(0.23 ±0.03)3 QgV^ 

ml = 0.8 GeV^ 
= 3.1 GeV 

mx = 3.87 GeV 

/v, = 0.405 GeV (37) 

The value of the angle a that defines the mixing be- 
tween the D°D*°, D°D*° and D+D*' D-D*+ has 
been obtained previously in Ref. [13, [13, HI] : 

a = 20° (38) 

For the mixing angle of two and four quark states, 9, 
we use the values that were obtained in the QCD sum 



rules analysis of the mass of the X and the decay mode 
X Jli)(n'n) [2|: 

61 = (9 ±4)°. (39) 

In the LHS of Eq. (^5]) . the unknown functions Gi{Q^) 
and Hi{Q^) have to be determined by matching both 
sides of the sum rule. In Fig. [2j we show the points 
obtained if we isolate the functions Gi{Q^) in Eq. 
and vary both and AP. The functions Gi{Q'^) [and 
consequently A{Q^), B{Q'^),C{Q'^)] should not depend 
on Af^, so we limit our fit to a region where the func- 
tion is clearly stable in to all values of Q^. We 
can see in Fig. [5] that the regions of stability in 
for Gi(Q2) is 7.0 GcV^ < AP < 8.5 GeV^ for G2iQ^) 
is 6.5 GeV^ < AP < 7.5 GeV^ and for GsiQ^) is 
8.0 GeV^ < M2 < 9.0 GeV^ 
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In Fig. ([2]) we show, through the dots, the QCDSR re- 
sults for the functions Gi{Q^) as a function of Q^. The 
form factors A{Q'^), B{Q^), and C{Q^) can be easily ob- 
tained by using Eas. ([5^ . and Since the cou- 
pling constants, appearing in Eq. ([M)) . are defined as the 
value of the form factors at the photon pole: = 0, to 
determine the couplings A, B and C we have to extrap- 
olate A{Q^), B{Q-^), and C{Q'^) to a region where the 
sum rules are no longer valid (since the QCDSR results 
are valid at the deep Euclidean region). To do that we 
fit the QCDSR results, shown in Fig. as exponential 
functions: 

G.(g')-5ie-^^'^'. (40) 

We do the fitting for sj^^ = 4.4 GeV and Uq^^ ~ 
3.6 GeV as the results do not depend much on this pa- 
rameters. The numerical values of the fitting parameters 
are shown in the Table ID 

From Fig. ([3]) we can see that the dependence of 





Gi 


Ga 


Ga 


52 


0.056 GeV=* 
0.25 GeV"^ 


-0.0069 GeV 
0.365 GeV"^ 


-0.013 GeV^ 
0.41 GeV"^ 



TABLE I. Results for the fitting parameters. 



the QCDSR results for the functions Gi{Q'^) are well re- 
produced by the chosen parametrization, in the interval 
2.5 GeV^ < < 4.5 GeV^ where the QCDSR are vahd. 

Using Eqs. ^ and (gOl) and varying 9 in 

the range 5° < 6* < 13° we get: 

A = A(Q2 = 0) = 18.65 ± 0.94 ; 
A + B = {A + B){Q'^ = 0) ^ -0.24 ± 0.11 ; 

C = C(Q^ = 0) = -0.843 ±0.008. (41) 
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The decay width is given in terms of these couphngs 
through [26} : 



result. 



{A + Bf + ^{A + C) 



a V 

r(x^j/V.7)=3^ 

(42) 

where p* — (m^ — m^)/(2mx) is the three-momentum 
of the photon in the X rest frame. To compare our re- 
sults with the experimental data shown in Eq. ([1} we use 
the result for the decay width of the channel J/^/itt+tt", 
obtained in the Ref. |24l |. which was computed in the 
same range of the mixing angle and with the same an- 
gle a = 20'': r{X -> J/?/- tttt) = 9.3 ± 6.9 MeV. We 
get 



TjX ^ J/ib 7) 

T{X J/tP TT+TT-) 



= 0.19±0.13. 



(43) 



which is in complete agreement with the experimental 



V. CONCLUSIONS 

We have presented a QCDSR analysis of the three- 
point function of the radiative decay of the X{3872) me- 
son by considering a mixed charmonium-molecular cur- 
rent. We find that the sum rules results in Eqs. are 
compatible with experimental data. These results were 
obtained by considering the mixing angles in Eq. (jlOp 
and © with the values a = 20° and 5° < 6^ < 13°. The 
present result is also compatible with previous analysis of 
the mass of the X state and the decays into J/'^tt'^tt+tt^ 
and J/i/jTT+TT" [iJl, since the values of the mixing an- 
gles used in both calculations are the same. It is im- 
portant to mention that there is no free parameter in 
the present analysis and, therefore, the result presented 
here strengthens the conclusion reached in Ref. [IJl that 
the X(3872) is probably a state with charmonium and 
molecular components. 
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